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Abstract: 

Large volcanic eruptions, in addition to the well-known effect of producing 
global cooling for a year or two, have been observed to produce shorter- 

term responses in the climate system involving non-linear dynamical proces- 
ses. In this paper, we use the ECHAM2 general circulation model forced with 

stratospheric aerosols to test some of these ideas. Run in a perpetual- 
January mode, with tropical stratospheric heating from the volcanic aero- 
sols typical of the 1982 El Chichdn eruption or the 1991 Pinatubo eruption, 
we find a dynamical response with an increased polar night jet in the Nor- 
thern Hemisphere (NH) and stronger zonal winds which extend down into the 
troposphere. The Azores High shifts northward with increased tropospheric 

westerlies at 60°N and increased easterlies at 30°N. Surface temperatures 

are higher both in northern Eurasia and North America, in agreement with 

observations for the NH winters or 1982-83 and 1991-92 as well as the win- 

ters following the other 10 largest volcanic eruptions since 1883. 
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(Lough and Fritts, 1987) and northern Europe (Groismann 1985, 1992). A 

reanalysis of global temperature data (Robock and Mao 1992, submitted) 

reveals winter warming over large parts of Eurasia after all major volcanic 

eruptions of the last century since the eruption of Krakatau in 1883. The 
question thus arises, whether there are processes other than the direct 

effect of radiation which determine the climate response to volcanic aero- 
sol during winter. One possibility is the dynamic coupling of stratosphere 
and troposphere, as has been shown with simple linear models (Geller and 
Alpert 1980, Schmitz and Grieger 1980). Numerical experiments (Boville 

1986) showed for instance reduced meridional heat transport in the tropo- 
sphere in the NCAR Gimate Community Model for conditions typical for stra- 
tospheric westerly circulation (cold polar and warm tropical stratosphere). 

We therefore conducted a GCM experiment in which, in addition to the reduc- 
tion of solar radiation at the top of the model atmosphere, heating rate 

anomalies were added in the stratosphere producing a rather realistic tem- 
perature anomaly field in the stratosphere. We want to study the atmosphe- 

ric response to the volcanic disturbance keeping the ocean fixed. Therefore 
any transitional effects from slowly varying parts of the climate system 

can be neglected. The atmospheric memory is at best two weeks. Since our 
main interest is in wintertime conditions, the model was integrated in a 
"Perpetual January"-experiment mode. Cumulative hydrological effects like 

the build-up of an enhanced snow cover and its effect on the Asian monsoon 
system are discussed in Graf (1992). 

Experimental Setup 

The stratospheric aerosol distribution observed in January 1983 in midlati- 

tudes of the northern hemisphere, ten months after the eruption of El Chi- 
chbn, is very similar to that seven months after Pinatubo in January 1992. 

The transmission of direct solar radiation lor wavelengths of X=0.3pm to 

2.8|im over Mauna Loa (sec Climate Diagnostics Bull., June 1992) is somewhat 
reduced in January 1992 (t=0.86) compared with January 1983 (x=0.89), and 
the column aerosol backscatter coefficient at ruby- wavelength (X=0.69 pm)at 

the same place is larger by a factor of two for January 1992. But, in midd- 
le latitudes of the northern hemisphere, lidar measurements of the aerosol 

column backscatter are comparable for both Januaries (Jager, pers. comm., 
see Figure 1). Later in 1992 the Pinatubo aerosol became more effective in 


4 




extinguishing solar radiation also in midlatitudes by a factor of 2 to 3 

compared to the El Chichbn aerosol. For our model simulation we used a 

prescribed aerosol distribution according to January 1992 observations, 

i.e. similar aerosol like in January 1983 for midlatitudes and stronger 

aerosol in low latitudes leading to a reduction in shortwave radiation as 

shown in Figure 2 (dashed line). Lads et al. (1992) found that the only 
important parameter for the net radiative forcing due to volcanic strato- 
spheric aerosol is the optical depth. Therefore, lacking further detailed 

information about the optical properties of the January 1992 aerosol (like 


SOLAR FORCING TOP O r ATMOSPHERE 


Figure 2: Reduction of solar radiation |W/m 2 | 
due 10 volcanic aerosol lor January 1 9 # 3 <so- 

HKii ,ne) . an<J 1992 (brHkcn lin d. The January 
1992 values are used as solar forcing in the 
experiments. 
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e.g. effective radii of the particles), we used model calculations based on 
El Chichbn to prescribe these. 

The climate model used in our computations is the ECHAM2 GCM, used at the 
Max-Planck-Institute fur Meteorologie for climate modelling research (Roe- 
ckner et al. 1992). We used the T21 version (horizontal resolution approxi- 
mately 5.6°) with 19 vertical levels, the full diurnal cycle and standard 

physics, including cloud-radiation interactions. Sea surface temperature 

(SST) 30(1 sea ice extent are prescribed from climatology. An investigation 

of the relative effects of anomalous SST assiociated with the 1991-92 El 

Nino and the volcanic aerosol is in preparation. 

Since the radiation scheme of the ECHAM2 climate model is unable to deal 
with variable aerosol concentrations, we decided to use the "anomaly forc- 
ing technique which has previously been used to study the climate response 

to the smoke from the burning Kuwait oil wells (Bakan et al. 1991). This 

technique allows the inclusion of any radiational effects of trace gases or 
aerosols in the model by the use of an additional radiation model, which 
exphcitely treats the species under investigation. In this case we used 

the 5- Eddington model, Bakan (1982). with the extensions for aerosol and 

ozone from Schult (1991). We ran this model parallel to the original ECHAM2 
radiation code once with and once without the prescribed aerosol with the 
atmospheric conditions from the climate model. The flux differences of the 
two 5-Eddington computations arcthen added to the results of the original 

code. After six months of running this complete interactive radiation and 

circulation code we obtained mean monthly anomalies of the radiative fluxes 

and heating rates. The effect of the aerosol longwave radiation on the 
surface radiation balance was more than one order of magnitude less than 
the total effect. It will therefore (also in accordance with Lacis et al. 
1992) be neglected in the computations described below. 

Since the fully interactive model increases the compulation time for ECHAM2 
by a factor of five, we avoided its explicit application in order to be 
able to run a series of experiments long enough for statistical investiga- 
tions. Thus, we used the externally . computed mean heating rate anomaly and 
the mean reduction of the shortwave radiation (Figures 2a and 3) of the 
fully interactive model for a scries of perpetual January simulations. The 
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model was initiated with data from one arbitrarily chosen January 1 of the 
control experiment and integrated for 30 days. Doing this 60 times (each 
new 1st January being initialized with the fields of the 30th January of 
die former integration) we produced a series of 60 January simulations both 
for the disturbed (Experiment) and the undisturbed (Control) case. 


Heating Rates Kelvin/Day 






Figure 4 gives an overview of the mean thermal and solar net radiation 

anomalies at earth surface as the climate model produced for three diffe- 
rent experiments for January conditions: Perm3 is our Pinatubo experiment, 
SST-2 is an El Nino experiment using the observed tropical (30°N to 30°S) 

sea surface temperature anomalies for January 1983, and Pcrm6 is the combi- 
ned El Nifio / Pinatubo experiment which will be discussed in a succeeding 
paper. Anomalies of surface longwave as well as shortwave radiation are in 
the same order of magnitude for all experiments. In all cases there is net 

increase of longwave radiation in middle and high latitudes. The pure El 
Nifio experiment only leads to a meridional shift of cloud areas and there- 

fore to a global mean zero effect of the shortwave anomaly, while the vol- 
canic aerosol reduces solar radiation in the order of 5 to 7 W/m 2 . 


NE THERMA. RADIATION SFC. ANOMALY NE~ SOLAR RADIATION S r C. ANOMALY 



F, ?r ^ et , anomal - v ( W/m'I ai the surface for thermal (Fig. 4a) and solar 
(l-ig. 4b) radiation for different Permanent January Experiments. PERM3 is 
the actual volcano experiment. SST2 is forced by tropical sea surface 

temperature anomalies as observed in January 1983. and PERM6 is a combi- 
nation of volcanic and SST forcings. 


We inspected carefully the results for possible trends in the zonal mean of 

temperatures and gcopotcntial and found that there were no strong trends 
detectable after few months. To check for the stability of any signal, we 

performed analyses with different subsets consisting of 30 Januaries each. 
So we chose blocks of 30 succeeding Januaries and scls of randomly extrac- 
ted Januaries, always excluding the first two months which show initial 

drift. All subsets gave the same principal patterns with only small varia- 
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lions in amplitudes and patiem. Here we will concentrate our discussion to 

the last half of the ensemble. The permanent January simulations are not 
independent from each other. This must be considered in the statistical 

tests by reducing the number of degrees of freedom. The effective number of 
degrees of freedom was computed for every gridpoint using the autocorrela- 
tion function. The resultant Bartels' number B (Bartels, 1935) has a global 
mean value of B = 1.8 (B = 1.4) lor the 2m temperature (500 hPa height). 

This means that the elfective degree of freedom is Ncir = N / B. This is 
consistent with the experience that the atmospheric memory is at most two 
weeks. Trying to stay at the conservative side of estimates, we made the 
assumption that every second experiment January is independent. Consequent- 
ly for a set of 30 Januaries we use an effective degree of freedom of 15. 
This leads to a slight underestimation of the significance of the anomaly 

amplitudes over most parts of the globe. Only over Antarctica the Bands’ 
number is larger, indicating some trend of the model in this region. 


Model Results 
1. Stratosphere 

Although the model stratosphere has only a coarse resolution, the obtained 
results arc comparable with observations and other model studies. 


The zonal mean cross section of the simulated temperature anomaly (Figure 
3b) is in good qualitative agreement with the corresponding obsen/ations 
Irom December 1991 (Figure 5).Positivc temperature anomalies arc simulated 
and observed in low latitudes, and negative anomalies in the polar night 
region of the stratosphere. Although the heating rate anomaly is positive 

over the North pole (see Figure 3a), the pole cools due to a dynamic ef- 
fect; because of the intensified polar stratospheric vortex the meridional 

(poleward) heat transport decreases. The meridional temperature gradient 

between the tropics and the winter pole is .strengthened in the lower stra- 
tosphere. In January 1992 a minor stratospheric warming occured increasing 
the temperature in high northern latitudes by 5 to 6 K. Although stratosph- 
eric warmings in principle can be simulated by ECHAM2. it is not a predict- 
able feature and is not reproduced by the model in this case. 
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Figure 5: Observed stratospheric Decanter 
temperature anomaly from the 1982-91 mean 
after the eruptions of El Chichbn (1982) 
and Pinatubo (1991), data from Labitzkc 
(radiosonde, solid lines) and Christy (MSU 
broken lines). ’ 


Our results are also consistent with the study of Rind ct al. (1992), who 

used the GISS model with an increased number of model layers (23 layers, 

top at 0.001 hPa). They investigated the effect of a homogenous aerosol 

layer of optical depth x=0. 15 on the zonal mean conditions in the middle 
atmosphere. 

The dynamic effects in the stratosphere are primarily due to changes in the 

meridional temperature distribution. The mean zonal wind anomaly (Figure 

3c) shows an enhanced polamighi jet in the stratosphere with anomalies up 
to 10 m/s. These westerly anomalies penetrate also downward into the tropo- 
sphere resulting in a poleward shift of the upper tropospheric westerly jet 

of the northern hemisphere. 

2. Troposphere 

The model simulates a significant reduction of the variance of the geopo- 

tential height of the 200. 500 and 850 hPa levels over the whole North 
Atlantic region. This holds also for the variances of the zonal wind field 
and temperature (the latter mainly in the lower troposphere) for the comp- 
lete set of the experiment (Figure 6) as well as for subsets. This is espe- 

cially clear for the gcopotential over the Azores high, over Labrador and 
Southern Greenland. No significant signal is detectable elsewhere applying 
the F-test. or the variance is increased in the experiment (Figure (i). 
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Geopotential Height 850 hPa gpm 



Temperature 850 hPa Celsius 
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U-Velocity 850 hPa m/s 
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Figure 6 : Changes in local variance (F-icst) at the 850 hPa level for 
Gcojxncmial height (.op), zonal velocity (middle), and temperature (hoc- 
tom) for the Permanent January Experiment. Shading indicates significant 
f tat ween Expcrimcni and Control. Values larger (smalicr) than 

the Com'ml. ’ C Va " anCC m lhc cx P crimcm « larger (smaller) than in 




A local T-test was used to lest for the significance of the amplitudes of 
the zonal wind, geopotential height and 2m-temperature anomalies in diffe- 

rent sets of 30 experiment months. Here we will focus on the results for 
the Januaries 31 to 60, i.e. on the second half of the experiment. Princi- 

pially, the anomaly patterns were similar in the other sets of each 30 
Januaries. 

in addition to the t-test, and in order to test for the difference of the 

distribution functions of the Control and the Experiment, a recurrency 

analysis (von Storch and Zwiers, 1988) was also applied to several sets of 

Januaries. T-test and recurrency test give comparable results, especially 

concerning the areas of statistically significant signals. In each set the 

pattern of anomalies of the geopotential heights is similar throughout the 

troposphere. 

The largest amplitudes of the geopotential heights in all tropospheric 

levels occur over the North Atlantic, where the pattern indicates a north- 
ward shift of the tropospheric circulation (Figure 7). The Azores high is 

displaced northward by about 15° and the Icelandic low shifted towards 

Greenland. This pattern is significant ((3=0.05) following the T-test and 

recurrent in all sets as well as are the negative tropical geopotential 

anomalies, mainly in the middle troposphere. The most prominent signal is 

the anticyclonic anomaly in the area of the northward shifted Azores high. 

Thus, the area of strongest and recurrent circulation anomalies is the 

North Atlantic. This holds for the whole troposphere. The structure of the 

anomaly patterns throughout the troposphere suggests a barotropic response. 

It is suggested that the interaction of the strengthened stratospheric 

polamight jet with tropospheric lee effects of the Rocky Mountains is 

responsible for the strongest response of the circulation to the aerosol 

forcing over the Nonh Atlantic. The primary effect may be a change in the 

stationary planetary wave pattern. 

The westerlies over the North Atlantic arc shifted northwards resulting in 

westerly wind anomalies around 60" N and easterly anomalies at 30° N (Figu- 
re 8). This pattern is again stable throughout the whole troposphere. It is 

also combined with a significant reduction of the variance of the zonal 

wind component as is obvious Irom the F-icst. Accordingly the Nonh Atlan- 
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tic centers of action shift northward. 


Geopotential Height 200 hPa fram 



Geopotential Height 850 hPa gpm 
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Fij:ure 7: Anomalies of Hk- 

levels. Shading indicates 
differences between Control and 
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Figure X: Anomalies of ihc zonal wind component at different tropospheric 
levels. Shading indicates statistically significant (error chance 5 %) 
differences between Control and Experiment for months 31 to 60 on the 
basis ol a local i-tcsi, degrees of freedom reduced to 15. 
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Figure 9: Anomalies of the air icmperaiurc in ihc lower troposphere. Sha- 
ding indicates statistically significant (error chance 5%) differences 

between Control and Experiment for months 31 to 60 on the basis of a 
local t-icst, degrees of freedom reduced to 15. 



In the tropics, the anomalies of the zonal wind component change sign from 
the lower to the upper troposphere. In both sets over the eastern equatori- 
al Pacific in the lower (upper) troposphere easterly (westerly) wind anoma- 
lies are found. This indicates a slightly strengthened Walker circulation 

over the eastern equatorial Pacific during winter. 

The meteorological parameter which normally is of interest in connection 
with violent volcanic eruptions is the surface air temperature. Figure 9 

(lower panel) shows the mean anomaly of this parameter, again for the last 
30 Januaries of our experiment. The main features, however, are similar 
also for other subsets of the 60 January integration after few months of 
initial drift. The T-test shows a major significant area of cooling over 

North Africa and the Middle East, strong cooling over Greenland and parts 

of China. In the southern hemisphere the only significant signals are a 

cold southern tip of South America and Antarctic cooling. Significant war- 

ming only occurs in the Barents Sea. Here one of the model s restrictions 
is apparent, i.e. the fixed sea ice coverage and temperature. This in part 

prevents the development of positive anomalies of the 2m temperature. In 
the free troposphere (e.g. in 850 hPa. Figure 9, middle panel), where the 

surface conditions are less important and the transport processes become 
dominant, the warming is significant over large parts of the Arctic ocean 
north of Siberia and the warm anomaly reaches further south over Asia. 

Tropospheric Observations 

In order to compare our simulation of the impact of the Pinatubo (or an 

other similar volcanic eruption) stratospheric aerosol on climate with 
global observations data, we used the 500 hPa geopotential height anomalies 
as analyzed by the US National Weather Service, the satellite obseived 

lower tropospheric temperature from the Microsounding Unit (MSU) (see Spen- 

cer and Christy. 1990). and (to a certain degree) synoptic observations 
collected by the Sccwciicramt of the Dcutscher Welterdicnsi. The last data 
set is based on a non-uniform distribution of stations with large gaps over 
the oceans. This gives serious problems when the data are interpolated to a 

coarse map such as the T21 grid. Thus these data are used only for qualita- 

tive comparisons. 

The MSU temperature data use the temperature dependence of the microwave 
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emission of oxygen molecules. For the tropospheric signal, calibrated MSU 

channel 2 (53.74 GHz) brightness temperatures are combined with a weight- 

ning function describing the vertical sensitivity to thermal emission by 

molecular oxygen in the atmosphere. The main pan of the signal we used is 

from the layer between the earth surface and 400 hPa (with a peak at 700 

hPa), i.e. from the lower pan of the troposphere rather than from the 

earth surface. The spatial resolution of these observations is in the order 

of 10°, even though the data are published on a 2.5° grid (Spencer et al , 

1990). 

In order to compare with previous violent volcanic eruptions, the surface 

temperature data set from Jones et al. (1988) was also used (see Robock and 

Mao, 1992 submitted). 

An inspection of the mean December 1991 through February 1992 500 hPa heig- 
ht anomaly (Figure 10) shows some similarity between our simulation and 
observed anomalies: The strongest anomalies occur in the Atlantic region. 


DJF 1991-92 



Figure 10: Anomaly of the 5(X) hPa gcopoicnlial height as analyzed by the US 
NWS. mean of December 1991 through February 1992. Anomalies calculated 
from mean 1982-1990. 
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The simulated anomalies are generally weaker than the observed ones and a 
tendency to a westward displacement is characteristic. Such zonal (west- 
ward) displacements (compared with observations) of the mean circulation 

features are common in low resolution GCMs, resulting mainly from the 

smooth orographic representation (Sausen 1992. pers. comm.). The negative 

anomaly over Greenland exists in both fields (for the simulated anomaly see 
Figure 6), while the observed center of the positive geopotential anomaly 

is shifted eastwards to England. Both, in the observations and in the simu- 
lation, a trough develops over Eastern Europe / Western Asia, which also 
has a more westerly position than in the observations. The pattern anomaly 

correlation (Table 1) is between 0.37 for the globe and 0.40 for the north- 
ern hemisphere if we use the mean pattern of the experiments 03-60. As 

stated above, subsets have shown the same patterns. This means that only 

about 15% of the total North hemispheric observed spatial variance of the 

500 hPa field of the middle troposphere in winter 1991-92 can be explained 
by the volcanic impact alone. Much better results we get for the western 
hemisphere (180°W to 0°) where this value can be improved to 49% (r=0.70). 

Table 1: Pattern correlation coefficients (area weighted, multiplied 

by 100) between observed and simulated meteorological fields. A 

1:1 correlation means that no spatial shift is allowed, "opti- 
mal” correlation allows for 2 gridpoints shift in all directions 

to get a best fit. All fields were interpolated to the T21 grid 

before being analyzed. Mean of the experiments 03-60 used. 


Parameter 

global 

N-hcmisph. 

90°W-90°E 

0°- 1 80”E 

0°-180°W 

MSU-Temp. 

1:1 

18 

57 

39 

54 

60 

optimal 

29 

57 

62 

56 

63 

500 hPa 
1:1 

35 

39 

1 

-29 

68 

optimal 

37 

40 

41 

-46 

70 

Siberia the 

patterns 

of 500 hPa 

geopotential 

height 

anomalies do 


agree well because of a simulated high over North Siberia which is not 
observed in winter 1991/92. 

For the lower tropospheric temperature (Figure 11), reasonable agreement 
between simulations and observations in winter 1991-92 is found for the 
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northern hemisphere again. The "optimal" pattern correlates are between 

t- 0.56 and t=0.63 for different subareas of the Northern hemisphere (Table 
1). resulting in x=0.57 for the whole northern hemisphere. This means that 
about 35 to 40% of the total spatial variance of the observed 2m temperatu- 
re can be explained alone by the volcanic forcing. The model reproduced 
especially well the cooling over Greenland and North Africa. Some of the 
warming of the middle and high latitudes of the northern continents, which 

is seen in both patterns, is also modelled at the right place, although not 

significant with respect to the amplitude because of the high variability 

of the temperature in these latitudes in the control run. 



The good agreement between the model simulation and the observations in 
winter 19991/92 still could be coincidental. This might be excluded if 

similar anomaly patterns evolved after other volcanic eruptions as well. 


An examination of the near surface temperature data set from Jones et al 
(1988 and 1991. updated. Jones, personal comm.) after the twelve most vio- 
lent volcanic eruptions of the last century (Figure 12) also shows good 

agreement with our simulations. For the six tropical volcanoes (Krakatau. 
Sta Maria. Agung, Fucgo, El Chichdn and Pinatubo) the anomalies are for the 
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uor i.ur race temperature Mnomuiy j 

After 12 Largest Volcanoes 1883-1992 





180 120W 60W 0 60E 120E 180 


Figure 12: Observed December through February mean temperature anomaly 

afier the twelve largest volcanic eruptions of the last century (Data from 
Jones, 1988. u|xl*aied, Jones, pers. comm.). 


following boreal winter, for the six high latitude volcanoes (Novarupta, 

Tarawera, Bandai, Ksudach, Cerro Azul and Bezymianny) anomalies are for the 
second winter after the eruption to allow time for the aerosols to reach 
the tropical stratosphere. The anomalies arc calculated from the 1951 to 
1980 mcanas differences from low pass (cut-off 10 years) filtered data. The 
pattern after every volcano is vinually the same as the mean effect shown 
in Figure 12, except some differences over North America for winters with 
and without El Nifios. The best pattern correlation between the composite 
signal of post eruption winter temperatures and our simulation is found for 
the sector between 90°W and 90°E on the Northern hemisphere (x=0.35). There 
is a zero pattern correlation for the western half of the North hemisphere 
and only x=0.21 for the eastern half. 


Conclusions 

The perpetual January GCM simulation has shown that the stratospheric aero- 
sol produced by a violent and sulfur-rich volcanic eruption can alter the 
stratospheric and tropospheric climate conditions significantly. Using the 
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“y f ° rCmS ' ' eCl,ni<, “ e forcing due to the aeresol 

‘ Ch » ««*- in ihc original OCM, can be simulated in acconTe 

sTl al 'can' beT l ^"’ UlilU ° nS a " d 0bSe,vaUons - Aspheric wi„ Br 

ta * in the . “ ' SI? "' nCara Warain S of ,hc ae ™ol coraaining 

layers ,n the tropics a„ d midlatitudes due „ abs0[p[io „ of radia[ive ^ 

gy is wanning increases the poleward meridional temperature gradient 
teads to an increased polar nigh, jet. reduced meridionai beat transpon 

therefore a cooling tendency in the stratosphere over the winter pole. 

in this study the tropospheric signal was found to be meridionally non- 

unifomv It ,s strongest over the Nonh Atlantic and consistsof enhanced 

westwmds over high latitudes (aiound 60°N>. strengthened zona, beat and 
moisture tntnspon ,n higher latitudes and a northward shift of the Azores 

' ' afeaS ° f slro "S cs ' an< ’ mali « a 'so show significantly reduced 

vanance of the parameters. 

A temperature anomaly field was found in the lower troposphere, which <i„ 

die main features) is very similar to obseivations during me winter ,99,. 

- imilar anomalies were observed also during other post volcano win- 

ters. There ,s a significant cooling over Greenland and Nonh 

nca/ Middle East. In agreement with observations some warming tendency 

Ith a maximum amplitude of about 2 K is seen in the mean results over 

orth Amenca and norihem Eurasia. This positive tempentture anomaly can 

not e proved to be statistically significant because of the large natural 

variability of surface temperature ,n these latitudes. This warn, anomaly is 

supponed by observations following historic eruptions (Robot* and Mao 

99.) and „ determines me high pattern coirelation found between Simula- 
uon and observations. 

The simulated and observed temperature anomalies following a violent 

VO came eruption ,n pan contradict the findings based on linear multiple 

ZT ,™' yS,S ' Cre “ and Sd, " nWi “ C l99l) ««<» balance models 

(Robock 1984,. The results underline the importance of nonlinear and 

dynamic responses of the atmosphere to external forcing. 
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